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Relativistic heavy ion collisions — Where are we now? Where do we go? 
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Various aspects of the current status of ultrarelativistic heavy ion colhsions are reviewed. 
Perspectives of heavy ion physics in the future are given as weU. 

1. INTRODUCTION 

The construction of RHIC (Relativistic Heavy Ion Collider) at BNL will be completed by 
CN I the time these proceedings are out and it will start accelerating gold nuclei at 100 GeV/A. 
^ I What is the purpose of ultrarelativistic heavy ion collisions? The reader may have dreams 
^Q • such as i) study of high density matter, ii) production of the quark-gluon plasma, iii) 
^ ■ quasi-reproduction of the big bang, iv) understanding of the history of the universe, and 
Q>^ . so forth. These are actually what will be pursued at RHIC, and so these are not dreams 
any more. We, however, have to turn our attention to the other side of reality as well: 
i) the lifetime of the system created in ultrarelativistic heavy ion collisions is very short, 
Q ' ii) the system is not static, iii) observables are dirty, in other words, the interpretation of 
^ ■ observables is not straightforward, iv) there are so many models that claim to describe 
the results of ultrarelativistic heavy ion collisions successfully, and so on. 

Thus, in order to fully understand ultrarelativistic heavy ion collisions and QCD, it is 
r> I really necessary for theorists to attack the following challenges, which were summarized 
S: by Matsui at Quark Matter 97 [|I1: 

(i) Compute, as best as one can, expected properties of dense matter and its phase 

structure and make predictions for signals of new states of matter. 

(ii) Interpret the data from current fixed-target experiments and identify signals of new 

physics, if any, from backgrounds of "old physics" . 

(iii) Describe formation and evolution of dense matter in nuclear collisions and estimate 

physical conditions to be achieved at future collider experiments. 

Since the first half of item (i) will be attacked by Yoshie 0, I will concentrate on the 
second half of item (i) and items (ii) and (iii). I will discuss mainly the recent results of 
heavy ion colhsions at CERN SPS at E^^, = 160 - 200 GeV/A. 

2. WHERE ARE WE? 

The system created in ultrarelativistic heavy ion collisions is not static. Even if 
hot/dense hadron matter or quark-gluon plasma is created, it cannot last for long time. 
The system immediately starts to expand and cool down. Even if the quark-gluon plasma 
phase is produced initially, it is soon converted back to the hadron matter, which freezes 



*This research was partly supported by Grant-in- Aid for Scientific Research no. 10740112 of the Japanese 
Ministry of Education, Science, and Cuhure. 



out in a time scale of a few tens of fm. A schematic diagram of the time evolution in 
ultrarelativistic heavy ion collisions is shown below. 
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Figure 1. Schematic time evolution expected in ultrarelativistic heavy ion collisions. 



This kind of explanation is often seen in the literature. However, the reader may wonder 
if this is indeed what is realized in ultrarelativistic heavy ion collisions. Without evidence, 
the whole picture is a mere conjecture. In the following, I will answer this question with 
recent experimental results and theoretical inference, and support this picture. 

2.1. Are hot systems ever produced in ultrarelativistic heavy ion collisions? 

Probably the reader is most interested in the possibility of the production of the quark- 
gluon plasma in ultrarelativistic heavy ion collisions. Nevertheless, I would like to begin 
with a more fundamental question; are there thermalization and collective motion in 
ultrarelativistic heavy ion collisions? It is because the phase transition to the quark- 
gluon plasma is a result of interaction. If no bulk interacting system is created, no phase 
transition takes place. It is logically possible that two nuclei collide on each other and go 
through each other without leaving a region with high energy density. However, recent 
experimental data tell us that approximately thermalized systems are indeed created in 
ultrarelativistic heavy ion collisions and that the systems go through collective expansion. 
The evidence includes the result of direct measurement of flow [^ , measurement of the 
interference of identical particles 0] , and so on. Here I will discuss the transverse mass, 
mrp, distribution of final state hadrons [^. 

The transverse mass, mj-, is defined by m^ = m^ +Pt, where m is the mass of the parti- 
cle and pt is its transverse momentum. The beam direction is defined as the longitudinal 
axis. The rriT distribution of final state hadron, i, at mid rapidity is well- approximated 
by the following form: 
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(1) 



Ti has been called temperature or slope parameter. It is known that Tj is a function of 
particle mass and well-approximated as 



Ti = a + brrii 



(2) 



where a and b are constants dependent on colliding nuclei, collision energy, and event 
class such as central or peripheral. This mass dependence of Tj has the following simple 



interpretation: a = Tf and b = (f ^), where Tf is the temperature of the system at freeze- 
out and V is the flow velocity at freeze-out. This relation is derived by assuming that 
all particles are locally in thermal equilibrium (not necessarily in chemical equilibrium) 
and collectively expanding at freeze-out and that the freeze-out time is independent of 
particle species. Thus, the experimental observation, Eq. (|^), is not inconsistent with the 
formation of hot thermalized system. 

Note that, as I emphasized above, this is merely one of the data that support the for- 
mation of hot systems. If this were the only evidence, the formation of hot systems cannot 
be concluded, since even in high energy pp collisions the rriT distribution is exponential 
[^. By combining various observables, the state of the system is deduced. This is the 
heavy ion way of inference. 

2.2. What kind of matter is created? 

We have learned that some kind of interacting matter is created, at least transiently, in 
ultrarelativistic heavy ion collisions. The next question is what kind of matter is created. 
One of the ways to get clues to it is to measure J/ip yield. 

As is well known, the suppression of J/i/j yield was originally proposed as a signature 
of the formation of the quark-gluon plasma 0]. The idea was that if the quark-gluon 
plasma is created, J/ip cannot form because the potential between a cc pair is Debye 
screened. However, this is not the only process that suppresses J/ip yield. Processes 
such as J/ip + N ^ DDN, can contribute to J/ip suppression. Thus, J/ip suppression is 
not necessarily a signature of the quark-gluon plasma by itself. Until PbPb experiments 
began at CERN, J/tp data in pA and AB collisions had been successfully explained by 
the above process, i.e., J/ip absorption by the nucleon. In PbPb collisions at ii^iab = 160 
GeV/A, however, it was found that as the transverse energy Et of final state hadrons 
increases, the J/ip yield, more strictly -B^^cr(J/'i/')/o"(Drell — Yan), drops suddenly at an 
Et 0. Hadronic scenarios failed to explain this behavior Q]. However, the nature of this 
sudden drop is not clear yet, although there are a lot of attempts to explain the behavior 
p!0| . In particular, I remark that even if the phase transition from the hadronic phase to 
the quark matter is of first order, it cannot lead to the sudden drop of the J/ip yield in 



such a naive way as discussed in ||rT|. It is because the energy density jumps at a first 
order phase transition. Thus, as Et is increased, part of the system begins to become 
the quark matter at E^i, and the portion of the quark matter increases gradually. When 
the whole volume has become the quark matter, E^ must have become much larger than 
Eti- As a result, the sudden drop in the J/ip yield does not correspond to the sudden 
formation of the quark-gluon plasma at a first order phase transition. 

2.3. How is the matter being excited? 

We have seen that some new form of matter appears to be created in ultrarelativistic 
heavy ion collisions. Then, how is the matter excited quantum mechanically? Dilep- 
tons are suitable probes to study this problem. Leptons do not interact strongly and 
can be considered almost penetrative. Observed dileptons carry the information of the 
early hot/dense stage as well as later stages, and some dileptons are produced even after 
freeze-out, for instance, by the decay of vector mesons. Accordingly, if vector mesons are 
modified, it is expected to be observed with dileptons, but not with hadrons. 

The first evidence of hadron modification was brought by CERES Collaboration at 



CERN SPS ]T^. They first measured dilepton and meson yields in pBe and pAu collisions, 



and found that the dilepton spectra in those reactions are explained solely by the decay 
of final state mesons. This means that no long-lived fireball is created in pBe or pAu 
collisions. They have, however, shown that the dilepton spectra in SU collisions at -Eiab = 
200 GeV/A cannot be explained only by the decay of final state mesons. Later, it was 
found that the data cannot be reproduced without taking account of hadron modification 
in the hot phase created in the collisions |ll3| , p |. 



Thus, hadrons are modified in medium. Two scenarios are often compared: mass 
shift and collisional broadening. The two scenarios are often considered two different 
scenarios. This is, however, not the case. First, I point out that the term, mass shift, is 
quite confusing. What is observed with dileptons is not the pole masses of vector mesons 
but the spectral function. In general, for the vector Heisenberg operator J^{x,t) the 
polarization tensor n^j,(go, (f) is defined by 

U^, = t fd'xe''^^{TJ^{x,t)Jl{6,0))T, (3) 



where (■ ■ ■)t indicates thermal average at T. For simplicity, let us set q = 0. The spectral 
function p{qo) is related to the polarization tensor. 

The physical significance of the spectral function is that the dilepton production rate at 
T, {dN^i/d^xd'^q)T, is related to the spectral function without approximation as follows 



0^ 77J^. TT^Pilo)- (5) 



d^xd^qj. 



oPqo 



This formula is exact, independent of in what phase the system is. There is some confusion 
on the meaning of dilepton spectra. Some authors argue that observed dilepton spectra 
are different from those theorists calculate because final state interactions modify theoret- 
ically calculated spectra. This statement stems from the misunderstanding that masses 
calculated by theorists are the masses of quasi-particles at T obtained by diagonalizing 
the Hamiltonian. If this were the case, the effect of final state interactions may change 
theoretical predictions substantially and should be taken into account. However, the for- 
mula, Eq. (^), is exact, and so no further interface between experiments and theories is 
needed except purely experimental corrections. 

Collisional broadening is a universal phenomenon; it appears wherever collisions take 
place |T^ . Generally, hadronic effective theories calculate this part of hadron modification. 
The other type of hadron modification, mass shift or global shift of peaks in the spectral 
function, is special to QCD. In the QCD sum rules, this is expressed as follows |T^. The 
operator product expansion side of the polarization function is given by 

On 



|rfW^*(TJ,(x,t)J^t(o,o))^ = -Cologig^l + E ^' (6) 



3Q27 ^ '^^ ' ' ^ ' '- ^ o,-. , . ^^^Q 



where Q^ = —uj^., Ci^s are condensates, and q = has been assumed as before. The 
condensates are related to the spectral function by the following exact sum rules: 

^%{^)s-ds = -^s^+' + i-irCn+u n>0, (7) 

n + 1 

where sq is the perturbative QCD threshold. At finite temperature or density, the conden- 
sates change because of partial chiral restoration... etc.. The change is reflected to global 
shift of peak structure in the spectral function at finite temperature or density through 
the exact sum rules, (0), in addition to the trivial collisional broadening [l^. Thus, the 



two scenarios for hadron modification are not exclusive with each other. In QCD both 
mechanisms are indeed at work and should be taken into account. In particular, the sum 
rules, (|^), should be satisfied by every effective model as well. 

2.4. Possibility of non-equilibrium states? 

Heavy ion reactions take place within a finite time. The typical time scale is of the 
same order as that of strong interaction and the process does not necessarily proceed 
adiabatically. Therefore, there is plenty of potential room for non-equilibrium phenomena. 
This non-equilibration is not simply limited to phase space, but also expected in chiral 
space. 

One of such possibilities is disoriented chiral condensate (DCC). Rajagopal and Wilczek 
proposed a mechanism to create DCC domains called 'quench mechanism' ||2^. I use, as 
a model Lagrangian, the linear sigma model defined by 

C = ^d,(jy,d^<f), - ^(0' - v')' + Ha, (8) 

where (pi = {a, tt) stands for a vector in internal space; Ha is an explicit chiral symmetry 
breaking term due to the finite quark masses. 

The original idea of the quench mechanism for DCC formation is summarized as follows. 
First, chiral symmetry is restored in the central region in particle or nuclear collisions. 
Then, the chiral fields are assumed to decouple quickly from the heat bath. The chiral 
fields are thus left around the origin in chiral space, i.e., 0^ ~ 0, while the effective 
potential has returned to its zero temperature form. As a result, the chiral fields are left 
at the top of the so-called 'Mexican hat' effective potential and according to its initial 
condition at each spatial point, the chiral field (p rolls down the slope of the effective 
potential in an arbitrary direction toward the bottom of the potential in chiral space. 
If in a certain spatial region the chiral field collectively rolls down in approximately the 
same direction and acquires an expectation value different from that in the vacuum, it 
will result in a DCC domain in coordinate space. This is schematically shown in Fig. 2. 

Actually, the above explanation tells us why we can expect non-equilibrium states in 
chiral space, but does not tell us what is the mechanism to create 'domains', i.e., large scale 
structure. It is the amplification of low momentum modes caused by the mode instability 
in the Mexican hat effective potential. The equation of motion for the momentum k 
component of the pion field vr* is given in the mean field approximation, 

^ = [A(.^-(0^))-fc^W, (9) 
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Figure 2. Configuration of tlie cliiral field 
(a) in the vacuum and (b) in a DCC state. 



Figure 3. Distribution probability of R 
defined by Eq. (|TUp. 



where the fluctuations of the chiral fields were neglected and (0^) is the average of 0^ = 
0*0*. In this approximation, modes with k"^ < v^ — (0^) grow exponentially, while high 
momentum modes do not. As a result, a large scale structure is expected to emerge. 

If a DCC domain is created and all final state pions are emitted from the single domain, 
it can be shown, using the remaining 0(3) invariance, that the distribution probability of 
a quantity R defined by 



R 



N: 



■TO 



Nno + K+ + K- 



(10) 



where Nt,^, Nt^+, and A''^- are the numbers of the final state ttq, tt~^, and n , respectively, 
takes the following form ||21||: 
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1 
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(11) 



Note that P{R) diverges at i? = 0, but that the average of R takes a finite value, 1/3. 
This probability distribution is quite in contrast to that expected for normal incoherent 
emission, P{R) ~ S{R — 1/3) (See Fig. 3). This behavior is valid only for low pt pions, 
because DCC formation is caused by amplification of low momentum pion modes. 

Eq. d^) tells us that it is necessary to get the system cooled fast enough to have 
(i) spontaneously symmetry broken effective potential and (ii) discordance between the 
minimum point of the spontaneously chirally broken effective potential and the position 
of the mean field in order to have large and discernible domains. 

In central ultrarelativistic heavy ion collisions, the initial fluctuation is expected to be 
large and the typical time scale for the expansion of the system is also large. Thus, the 
system is likely to maintain equilibration at least in chiral space. The above conjecture 



has been also confirmed numerically by our group [^], assuming the Bjorken scaling in the 
longitudinal direction ||23|. This leads to the following conclusion: central ultrarelativistic 
heavy ion collisions are not the best place to look for DCC formation, contrary to the 
general expectation. It is not-so-central ultrarelativistic heavy ion collisions that should 
be collected to look for DCC formation. I also note that in not-so-central ultrarelativistic 
heavy ion collisions the axial anomaly is expected to bring about DCC domain formation 
aligned in real space, i.e., one above the reaction place and one below it, but misaligned 



in chiral space p4|. This conjecture can be tested by utilizing the technique to determine 
the reaction plane, which has been developed in flow analysis. 

At the moment, only one heavy ion experiment has reported the result of a DCC search 
[p5| . Their result was negative. However, they used only central collisions and did not 
imposed a cut on px- As I discussed above, it is necessary to select non-central but non- 
peripheral collisions and low pt pions in DCC hunts. Until experiments with such cuts 
are carried out, the possibility of DCC formation in ultrarelativistic heavy ion collisions 
will remain an open question. 



3. WHERE DO WE GO? 

RHIC is coming late this year (1999) at v^ = 200 GeV/A (A = ^^^Au). Hopefully, 
LHC will come early next century. While physical observables and methodology will 
remain essentially the same, energy will be much higher, momentum and invariant mass 
resolution will be much better (the e~^e~ invariant mass resolution is about 20 - 30 times 
better at 1 GeV than the current CERES setup at CERN SPS), and acceptance will be 
much larger. The quark-gluon plasma, whose formation is not decisive at the moment, will 
be created with much larger probability. Since the decrease of the temperature is slowed 
near a phase transition irrespective of the order of the transition, if the phase transition 
occurs, it will be possible to separate dileptons from the (almost) constant temperature 
period near the critical temperature as a secondary peak or shoulder structure in the 
invariant mass distribution [^ . 

The system created in ultrarelativistic heavy ion collisions is very complex. In order 
to obtain profound understanding of the small short-lived non-static and potentially non- 
equilibrated system, it is mandatory to compare as many observables as possible and find 
correlations among them. The comparison should not be limited to event classes at a 
fixed collision energy, but should be among events at different collision energies. So far, 
there has been a large gap in collision energies of heavy ion experiments between i^iab = ~ 
12 GeV/ A (BNL AGS) and -160 GeV/A (CERN SPS). The gap between CERN SPS and 
RHIC is considerably large as well. From this viewpoint, the approval and inauguration 
of the JHF project, which is capable of accelerating heavy ions at E'lab = 50 GeV/A and 
lower, are eagerly awaited. It is also desirable to run RHIC at lower energies. 
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Figure 4. Energies of past, present, and future high energy heavy ion machines. 



Finally, I make comments on event generators. Quite a number of event generators exist 
in the market |2^ . Most of them are classical, although Pauli blocking is generally taken 



into account. This is acceptable as a first step. However, it is not acceptable that most 
models do not include quarks or gluons and that many models do not include secondary 
collisions. This is completely contrary to the general expectation that RHIC physics is 
described by quark and gluon degrees of freedom and the fact that phase transition takes 
place owing to interaction among quarks and gluons. In this respect, perhaps the most 
ambitious attempt was Parton Cascade (or VNI) |^ by Klaus Kinder Geiger, although 
it is not without serious conceptual defects yet. He tragically perished in the Swiss Air 
Flight 111 crash on September 2, 1998, but I hope that this is not the end. For better 
understanding of the dynamical features of RHIC and LHC physics, the direction initiated 
by him should be taken over and further pursued. 

REFERENCES 

1. T. Matsui, Nucl. Phys. A638 (1998) 19c. 

2. T. Yoshie, these proceedings. 

3. For example, J. Barrette et al., E877 Collaboration, Phys. Rev. Lett. 73 (1994) 2532. 

4. I. G. Bearden et al., NA44 Collaboration, Phys. Rev. Lett. 78 (1997) 2080. 

5. H. Appelshauser et al., NA49 Collaboration, Nucl. Phys. A638 (1998) 91c. 

6. U. Heinz, Nucl. Phys. A638 (1998) 357c. 

7. T. Matsui and H. Satz, Phys. Lett. B178 (1986) 416. 

8. M. C. Abreu et al, NA50 Collaboration, Phys. Lett. B410 (1997) 337. 

9. R. Vogt, Phys. Lett. B430 (1998) 15. 

10. K. Geiger and B. Miiller, Heavy Ion Phys. 7 (1998) 207; E. V. Shuryak, Nucl. Phys. 
A638 (1998) 207c. 

11. D. Kharzeev, M. Nardi, and H. Satz, |hep-ph/ 9707308 . 



12. G. Agakichiev et al., CERES Collaboration, Phys. Rev. Lett. 75 (1995) 1272. 

13. C. M. Hung and E. V. Shuryak, Phys. Rev. Lett. 75 (1995) 4003. 

14. G. Q. Li, C. M. Ko, and G. E. Brown, Phys. Rev. Lett. 75 (1995) 4007. 

15. H. A. Weldon, Phys. Rev. D 42 (1990) 2384. 

16. C. Gale and J. Kapusta, Nucl. Phys. B357 (1991) 65. 

17. For example, see, A. Gallagher, Phys. Rev. A 12 (1975) 133. 

18. T. Hatsuda, Y. Koike, and S. H. Lee, Nucl. Phys. B394 (1993) 221. 

19. M. Asakawa and C. M. Ko, Phys. Rev. C 48 (1993) R526. 

20. K. Rajagopal and F. Wilczek, Nucl. Phys. B404 (1993) 577. 

21. J. -P. Blaizot and A. Krzywicki, Phys. Rev. D 46 (1992) 246. 

22. M. Asakawa, Z. Huang, and X. N. Wang, Phys. Rev. Lett. 74 (1995) 3126. 

23. J. D. Bjorken, Phys. Rev. D 27 (1983) 140. 

24. M. Asakawa, H. Minakata, and B. Miiller, Phys. Rev. D 58 (1998) 094011. 

25. M. M. Aggarwal et al., WA98 Collaboration, Phys. Lett. B420 (1998) 169. 

26. M. Asakawa and C. M. Ko, Phys. Lett. B322 (1994) 33; Phys. Rev. C 50 (1994) 3064. 

27. For example, see, Y. Pang, Nucl. Phys. A638 (1998) 219c. 

28. K. Geiger and B. Miiller, Nucl. Phys. B369 (1992) 600. 



